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ABSTRACT 

We analyze the constraints imposed by experimental upper limits on electric dipole mo- 
ments (EDMs) within the Maximally CP- and Minimally Flavour- Violating (MCPMFV) 
version of the MSSM. Since the MCPMFV scenario has 6 non-standard CP-violating 
phases, in addition to the CP-odd QCD vacuum phase ^qcd, cancellations may occur 
among the CP-violating contributions to the three measured EDMs, those of 205 T1, the 
neutron and 199 Hg, leaving open the possibility of relatively large values of the other CP- 
violating observables. We develop a novel geometric method that uses the small-phase 
approximation as a starting point, takes the existing EDM constraints into account, and 
enables us to find maximal values of other CP-violating observables, such as the EDMs of 
the Deuteron and muon, the CP-violating asymmetry in b — > S7 decay, and the B s mix- 
ing phase. We apply this geometric method to provide upper limits on these observables 
within specific benchmark supersymmetric scenarios, including extensions that allow for a 
non-zero 6> QC d- 
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1 Introduction 



The non-observation of electric dipole moments (EDMs) for Thallium ( 205 T1) [T], the neu- 
tron (n) [2], and Mercury ( 199 Hg) [3111] provide very tight bounds on possible new sources 
of CP violation beyond the Cabibbo-Kobayashi-Maskawa (CKM) phase of the Standard 
Model (SM). Prima facie, these bounds suggest that any CP- violating phases associated 
with new physics at the TeV scale are very small, posing a challenge to scenarios of new 
TeV-scale physics, such as supersymmetry (SUSY) [5], that contain many potential new 
sources of CP violation. The EDM challenge is compounded by the excellent agreement 
of present experiments with the CKM model [6], providing important complementary con- 
straints on the flavour structure of any new TeV-scale physics, as well as on its role in 
CP violation. On the other hand, the baryon asymmetry in the Universe (BAU) could be 
explained by TeV-scale physics, if it has substantial CP violation and realizes a sufficiently 
strong first-order electroweak phase transition in the early Universe [7]. 

This tension between a TeV-scale origin of the BAU and EDMs will be explored both 
directly and indirectly in the LHC era. Experiments at the LHC, notably LHCb [8], will 
soon be giving direct new information about flavour and CP violation at the TeV scale. In 
parallel to these direct explorations at the LHC, a new generation of precision low-energy 
experiments will play an important indirect role. These new precision experiments will 
place much stronger indirect constraints on the possible CP and flavour structure of models 
of TeV-scale physics. New experiments on the neutron EDM are underway and, if the 
proposed experiment searching for a Deuteron ( 2 H + ) EDM achieves its design sensitivity [HI 
[TO] , it will improve the existing bounds on possible CP- violating chromoelectric operators 
by orders of magnitude [TT) . 

In this paper we introduce a new geometric approach to the incorporation of EDM 
constraints on CP-violating models, showing how the maximal values of unmeasured CP- 
violating observables may be estimated in a systematic and reliable way. We illustrate this 
approach in the context of SUSY, regarding it as an archetype of the TeV-scale models 
that are (potentially) embarrassed by experimental constraints on flavour and CP violation. 
However, our geometric approach could easily be applied to other models, and indeed to 
other classes of observables besides those that violate CP. 

For illustration, we work within the minimal supersymmetric extension of the Stan- 
dard Model (MSSM), with SUSY broken softly at the TeV scale. We assume a generic 
version of this model with minimal flavour violation (MFV), driven by the fermion Yukawa 
couplings. As discussed in [TJ] and elaborated further in [131 E], this model has a total 
of 19 parameters, of which 6 violate CP. In the convention where both the superpotential 
Higgs-mixing parameter /i and the respective soft SUSY-breaking Higgs-mass term Bfj, are 
real, these are the three phases of the soft SUSY-breaking gaugino masses Mi^.3 and the 
three phases of the trilinear SUSY-breaking parameters A^u,e- This model was called in [TJ] 
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the Maximally CP- and Minimally Flavour- Violating (MCPMFV) scenario of the MSSM, 
or the MCPMFV SUSY model in short. In addition to the 6 MCPMFV phases, one may 
consider the CP-violating QCD vacuum phase 6*qcd- The specific question we study in this 
paper is how to find the maximum value of some other CP-violating observable, such as 
the CP- violating asymmetry in b — > 57 decay, Aqp, the phase in B s mixing, <pB a , or some 
other EDM that has not yet been measured (accurately), such as those for the Deuteron 
and muon, while implementing the available EDM constraints. 

In order to demonstrate the principle of our geometric approach, we first consider 
a toy warm-up problem involving a single constraint in a theory with three parameters. 
We characterize the subspace of parameters satisfying this constraint in the linear approx- 
imation, and then show how to identify the direction in this subspace that maximizes 
any given observable O. The generalization to multiple constraints in higher-dimensional 
spaces follows similar geometric principles, that may be described in the language of exte- 
rior products of differential forms [15] . This geometric optimization method of the so-called 
"cancellation mechanism" [16] may sound complicated, but its numerical implementation 
is straightforward. We emphasize that our geometric approach differs in principle from the 
naive scan method that is usually adopted in the literature [T7]. The geometric method 
proposed here provides an accurate parametric determination of the optimal cancellation 
regions where any given physical observable is maximized in the linear approximation. 
Hence, our geometric approach is exact, efficient and less computationally-intensive than a 
naive scan of a mult i- dimensional space. 

In the application of this approach to the MCPMFV SUSY model, we first select some 
benchmark points in the CP-conserving restriction of the model. We then evaluate the 
dependences of the relevant constraints (the EDMs of 205 T1, the neutron and 199 Hg) on the 
six CP-violating phases of the full MCPMFV model in the linear approximation, as well as 
the linear dependences of the CP- violating observables of interest (the Deuteron and muon 
EDMs, Aqp and 4>b s )- We then demonstrate numerically how each of these observables 
may be maximized in the linear approximation, taking into account the existing EDM 
constraints. We note that, at any specific benchmark point, the values of all CP- violating 
observables are bounded in magnitude, since the ranges of the CP-violating phases fa are 
all compact: fa G [0,27r). We also note that our approach is only approximate for large 
CP-violating phases. Nevertheless, the linear expansions give good estimates of the true 
maximal values of the CP-violating observables. 

We find that the linear approximations to the EDMs and other CP-violating ob- 
servables in the neighbourhoods of the MCPMFV benchmark points we study are quite 
accurate for CP- violating angles with magnitudes up to several tens of degrees. We confirm 
that the EDM-free directions in parameter space constructed using our geometric construc- 
tion yield values of the other CP-violating observables that are larger than those possible 
(in the linear approximation) along other directions in the space of CP-violating phases. 
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Along the optimal EDM-free directions, we find that values of the Deuteron EDM an or- 
der of magnitude larger than the prospective experimental sensitivity may be attained for 
acceptable values of the MCPMFV phases, and almost an order of magnitude larger still if 
the optimal geometric construction is extended to include the CP-violating QCD vacuum 
phase Oqcd- On the other hand, we find maximal values of the muon EDM that are below 
the likely experimental sensitivity in both the scenarios with and without the QCD phase. 
In the case of Aqp, we find values as large as 2% in the large tan/3 scenario. Given the 
constraint from B(b — > 57), however, Aqp cannot exceed the 0.1% level and so remains too 
small to be observed. Finally, the S s -meson mixing phase <pB s turns out to be close to the 
small SM value in both the scenarios studied. 

The layout of the paper is as follows. Section 2 presents our novel geometric approach, 
starting with the toy three-dimensional exercise and continuing to the six-dimensional case 
of the MCPMFV SUSY model. The implementation of this approach for a selection of 
benchmark points is described in Section 3, and our numerical results for Aqp, <Pb s and 
the Deuteron and muon EDMs are presented in Section 4. Section 5 extends the previous 
discussion to the seven-dimensional case including #qcd- Finally, Section 6 presents our 
conclusions and some suggestions for future work. 

2 Optimal EDM-Free Directions 

The current experimental upper bounds on the Thallium, neutron and Mercury EDMs 
put very strict constraints on the CP-violating parameters of the theory, such as CP-odd 
phases. However, they do not constrain all possible combinations of the CP- violating phases 
in models with many such phases, such as the MSSM. It is therefore important to develop a 
powerful approach for finding the optimal choice of CP-odd phases which maximize the size 
of a given CP- violating observable O, while remaining compatible with the present EDM 
constraints. Examples of such CP-violating observables for which the present experimental 
sensitivities are likely soon to be improved significantly include the Deuteron and muon 
EDMs, the CP asymmetry in b — > sj, and the phase in B s mixing. 

In this paper, we propose a geometric approach to the maximization of such observ- 
ables, which may also be applied to other analogous problems. For illustration, we first 
consider a simple three-dimensional (3D) example for a theory which has only three physi- 
cal CP-odd phases. Working in the 3D vector space defined by these three CP-odd phases, 
we show how to construct geometrically optimal 'EDM-free' directions in the small-phase 
approximation. We then generalize this approach to theories with more than three CP-odd 
phases, such as the MCPMFV SUSY model with its six new CP-odd phases, optionally 
including the QCD vacuum phase ^qcd. Using simple ideas from the calculus of differential 
forms [15] for such a higher- dimensional setup, we are able to derive the optimal EDM-free 
directions for any given CP-violating observable O that we wish to study. 
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2.1 A Simple 3D Example 



For the sake of geometric familiarity, we first consider the simple 3D example of a the- 
ory with just three physical CP-odd phases, represented by the 3D phase vector $ = 
($!, $2, ^3) • For illustration, we assume that we have a single very stringent EDM con- 
straint, which we denote by E. In the region of small phases, e.g., for |$| < tt/6, we seek to 
maximize the value of some specific CP-violating observable O, under the condition E = 0. 
Both the CP- violating observable O and the EDM constraint E are functions of i.e., 
O = 0(<&) and E = E(<&), and vanish in the limit of vanishing CP-odd phases: O(0) = 0, 
E(0) = 0- 

We make Taylor expansions of O and E in terms of the small phases, and keep only 
the linear terms in these expansions: 

O = $ • O , E = $E , (2.1) 

where we have defined O = VO, E = V-E 1 , and V = (d/d<f>i, d/d(p2, d/dcftz). The condition 
E — requires that the phase vector $ should lie in the plane orthogonal to E. In order 
to maximize the value of O, the phase vector $ should lie along the intersection of the 
plane spanned by the vectors O and E with this plane perpendicular to E, as represented 
schematically in Fig. [TJ Up to an overall sign, the solution is unique and is given by the 
double cross-product: 

= Ex(OxE). (2.2) 

Evidently, the condition E = is satisfied by construction in the small-phase approxi- 
mation, i.e., E = • E = 0. In this way, we can construct unambiguously the optimal 
EDM-free direction <fr*. 

It is straightforward to obtain the maximum value of the CP- violating observable O. 
In the small-phase approximation, this is simply given by 

O = <f)* <&* O = ±<p* tJ\0\ 2 - (O-E) 2 , (2.3) 

where 0* = and the carets on the vectors indicate unit-norm vectors. In the small- 
phase approximation, the largest possible value for the CP-violating observable O is ob- 
tained when E is perpendicular to O, in which case we have 

O max = ± <f>* |0| . (2.4) 

On the other extreme, if E happens to be parallel to O, then the CP-violating observable O 
vanishes, i.e., O = 0, in the limit that E = 0. Note that there is a twofold degeneracy in 

*Our approach can easily be extended to cases where some of the phases $1,2,3 approach another CP- 
conserving point different from zero, e.g., $1.2,3 — > ±tt, in which case the CP-odd phase vector $ represents 
the difference from this CP-conserving point. 
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Figure 1: Geometric construction of the optimal EDM-free direction in the small-phase 
approximation for a CP-violating observable O, subject to the EDM constraint E = 0, in a 
theory with three physical CP-odd phases $1,2,3- The optimal EDM-free direction is given by 
the CP-odd phase vector which is the intersection of the indicated plane perpendicular 
to the vector E with the plane defined by the vectors O and E. 

the optimal value of the CP-violating observable O, i.e., our geometric construction leaves 
O undetermined up to an overall sign. However, this twofold degeneracy is a consequence 
of the linear small-phase approximation, where quadratic and higher-order derivative terms 
in a Taylor series expansion, e.g., ViVjE, were neglected. These terms break this twofold 
degeneracy in general, and we include them all in our numerical analysis. 

2.2 A 6D Example: the MCPMFV SUSY Model 

The above geometric construction for the 3D example can be generalized to theories with 
more than three CP-odd phases, such as the MCPMFV SUSY model, which has six phasesQ, 
and including more than one EDM constraint. In the case of such a theory with six CP- 
odd phases, <fr is now a 6D phase vector, subject to three EDM constraints denoted by 
E a,b ' c = 0, corresponding to the non-observation of the Thallium, neutron and Mercury 
EDMs. As before, the task is to maximize a given CP- violating observable O, now satisfying 
simultaneously the three conditions: E a = E b = E c = 0. 

The generalization of the differential operator V to 6 dimensions is obvious, and it 
may be used to obtain in the small-phase approximation the four 6D vectors: E a ' b ' c = 
S7E a,b,c and O = VO. For simplicity, we assume that the four vectors E a,fe,c and O are 

'''We extend this discussion later to include the QCD phase #qcd- 
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linearly independent, i.e., there are no degeneracies between the EDM constraints and the 
observable O: we return to this issue towards the end of this Section. 



Given the above assumptions, the analogue of the single vector E in the 3D example 
is the triple exterior product 

A a(Sy = E[ a EpE^ , (2.5) 

where the Greek indices label the components of the vectors in the 6D space, i.e., a, (3, 7 = 
1,2, . . . , 6. The square brackets on the RHS of (12.51) indicate that the tensor A a/3l is 
obtained by fully antisymmetrizing the vectors E%, Ep and E'i in the indices a, /?, 7, i.e., 
A a /3-y = —A/3 ai = —A ai/ 3 etc. Borrowing a term from the calculus of differential forms, 
A a p 7 is a 3-form. 

Correspondingly, the analogue of the direction O x E, which determines the normal 
to the plane defined by O and E in the 3D example described above, is the 2-form 

B^u = £ pvXpoT O \ E a p E b a E^ , (2.6) 

where summation over repeated indices is implied and e^xpar is the usual Levi-Civita ten- 
sor generalized to 6D. In the language of differential forms, B pu is, up to an irrelevant overall 
factor, the Hodge-dual product between the 1-form 0\, representing the CP-violating ob- 
servable, and the 3-form A a ^. 

The components $* of the optimal EDM-free direction maximizing O can now be 
obtained from the Hodge-dual product of the 3-form Ap^s and the 2-form B^ v . Explicitly, 

E-y Bs 0\ Ep E a E C T , (2.7) 

where we have included an unknown overall normalization factor M . By construction, the 
6D phase vector <fr* is orthogonal to the three vectors E a,6,c , thus satisfying the desired 
EDM constraints, E a = E b = E c = 0, in the small-phase approximation. We observe 
that the magnitude <fi* = and hence the overall normalization factor A/", can only be 

determined by a numerical analysis of the actual experimental limits on the three EDMs. 
As in the 3D example, the maximum allowed value of the CP-violating observable O is 
given in the small-phase approximation by 

O = <f>* C K = ±AT \e^ aM& e^xpor O a O x E% E% E c 5 E a p E b a E c T \, (2.8) 

where the caret denotes the components of a unit-norm vector. As in the 3D case, quadratic 
and higher-order derivative terms with respect to the CP-odd phases will generically prefer 
a particular sign for the optimal value of O 0. 

As a consistency check of our geometric construction, one may verify that in the 
small-angle approximation the largest possible value, O max = ±0*|O|, is obtained when all 



■"■In some circumstances, the small-phase approximation may even break down in such a way that the 
maximal value of O is obtained for a value of <fi* less than the maximal value allowed by E. 
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the four 6D vectors, O and E a ' b ' c , are mutually orthogonal, exactly as in the 3D case. For 
example, if the vectors E a ' b ' c have non-zero components only in the 4, 5 and 6 directions, 
respectively, and the CP- violating observable O is a reduced 3- vector living in the 1, 2 and 
3 coordinates, it is not difficult to check that the expression (12 .8p gives the largest possible 
value O max , as expected. 

The above geometric construction can be extended to include further EDM or other 
strict CP-violating constraints on the theory. Specifically, in the MCPMFV SUSY model, 
we can afford to have at least two more linearly-independent constraints coming from 
different experiments and still be able to potentially get one large combination of CP-odd 
phases. If E d,e = are these two extra constraints, the maximum allowed value for the 
CP-violating observable O will then be 



O = ± Af e^afsys E a u E b a Ep E d Eg 



(2.9) 



where the normalization Af depends on the actual size of <fi*. Equation (12. 9p is nothing else 
than the simple generalization of the well-known triple-product to 6D. 

We can also allow for the possible presence of a non-zero strong CP phase #qcd i n the 
theory, in which case the corresponding CP-odd phase vector <fr becomes seven-dimensional 
(7D) in the MCPMFV SUSY model. In this case, four additional linearly-independent 
EDMs or other strict CP-violating constraints would be needed, in addition to the three 
limits from the present EDM experiments, in order to span fully the CP-odd space of the 
MCPMFV model, and so constrain the norm of the 7D CP-odd phase vector <fr as discussed 
in Section 5. 

In our geometric construction, an important role is played by the degree of degeneracy, 
or alignment, between pairs of observables, e.g., between O and the E a . For this purpose, 
it is interesting to know the cosine C*o,E a of the relative angle between their corresponding 
vectors O and E a , i.e., 

O E a 

Co * a = \ofW\ ' (2 ' 10) 

If C*o,E a = ±1, we have perfect alignment of the observables O and E a . In such a case, 
if E a vanishes, then so does O. On the other hand, if Co,E a = 0, the two observables O 
and E a are orthogonal, and hence can vary independently of each other. Using (I2.10p . one 
can obtain an upper bound for the optimal value of the CP-violating observable O in the 
small-phase approximation: 



101 < 



Max ( C? 



0,E a 



1/2 , 

|0|, (2.11 



where the index x labels the vectors E a,6 ' Cv " related to the different EDM constraints. 
Notice that the inequality (12. lip becomes an exact equality in the 3D example which we 
discussed before in Section |2jl [cf. (12. 3p . 
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Finally, after completing the description of this geometric approach, we note that it 
could be applied to many other phenomenological problems where there are M constraints 
on a theory with N > M parameters, and one wishes to determine the maximum value of 
some other observable quantity. 



3 Numerical Illustrations 

In this Section we illustrate the geometric approach introduced in Section [2]by constructing 
optimal EDM-free directions for some specific benchmark scenarios. For this purpose, we 
consider CP-violating variants of a typical CMSSM scenario with 

|M 1)2i3 | = 250 GeV, 

M 2 Hu = M 2 Hd =Ml = Ml = Ml = M\ = Ml = (100 GeV) 2 , 

|AJ = \A d \ = \A e \ = 100 GeV, (3.1) 

at the GUT scale, introducing non-zero CP- violating phases and varying tan /3 (Msusy)- We 
adopt the convention that $ M = 0°, and we vary independently the following six MCPMFV 
phases at the GUT scale: $i, $ 2? $3> $a u) $A d , and $A e - We note that the $1,2,3, and 
the fi parameter, $ M , are unchanged by the RG evolution at the one- loop level, whereas 
the phases of the trilinear couplings A u ,d, e at low scales could be significantly different 
from the values specified at the GUT scale. This scenario becomes the SPSla point [18] 
when tan/3 = 10, $1,2,3 = 0° and §A u ,A d ,A e = 180°. Our calculations of the EDMs include 
the two-loop diagrams mediated by the / y-H ± -W T and , ~f-W ± -W T couplings, which are 
summarized in the Appendix. Unlike [19], we incorporate the $3 dependence induced by 
gluino exchange in the one-loop H ± -u-d coupling, which becomes relevant in the region of 
large tan/3. 

In order to analyze this scenario, we first make Taylor expansions of the following five 
EDMs and two CP-violating observables: 

(fan/ 5 d>n/ dn XP > ^H g / ^Ig P j ^d/ ^d XP j dp/ > (3-2) 
AcAb -> sj)[%} , <p Bs = Arg (<4W =2 |£ s °>susy) [°] , 

where we choose the following normalization factors 

d*f p = 9 x 10- 25 e cm , ^ xp = 3 x 10" 26 e cm , d|f = 3.1 x 10" 29 e cm , 
d* xp = 3 x 10- 27 e cm , df p = 1 x 10~ 24 e cm. (3.3) 

In the cases of the EDMs of Thallium [T], the neutron [2], and Mercury [3111], we use the 
current experimental limits for normalization, and for the EDMs of Deuteron and muon 
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Figure 2: TTie quadratic fit to the Thallium EDM that is used to obtain the 6D vector 
E dT1 = V(c?Ti/^fi CP ) ^ n an expansion around (p a = 0° for the scenario ( fff. 1\) with tan /3 = 40. 



we use the sensitivities projected in [91(20] . The observables A CP and (p Bg are measured in 
percent and degrees, respectively. 

In order to obtain the vectors representing the EDM constraints and the observables 
in the six- dimensional CP-phase space, namely E = VE and O = VO where 

V Q = (d/d^d/d^d/d^d/d^d/d^A^d/d®^) , (3.4) 

we calculate the EDMs and the observables in ranges A$ = ±10° around each chosen CP- 
conserving point (p a , varying independently each of the six CP phases, and performing a 
quadratic fit in each case. As an example, Fig. [2] shows a quadratic fit to the Thallium EDM 
for the scenario under consideration taking tan f3 = 40. We observe that the CP-violating 
phase dependence is in fact linear in the region considered, to a very good approximation. 

In Figs. [3] and HI we show the absolute values of the components of the six 6D vectors 

E d ^ = V(d Tl /d*? p ), E^ = V(djd* xp ), E dH * = V(d Hg /4f); 

O dD = V(V^ XP ), O^V(Vdf p ), 0^=V(A CP /%), (3.5) 

around the CP-conserving points (p a = 0° and (p a = 180°, respectively, varying tan/3. 
We do not show 0^ Bs , since we find that this observable is too small to be detectable in 
the class of scenarios under consideration. The solid lines are for the components of the 
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Figure 3: The absolute values of the components of the three 6D vectors representing the 
present EDM constraints (upper), and those of the three 6D vectors representing other CP- 
violating observables (lower) in expansions around the CP-conserving point ip a = 0° as 
functions of tan (3 for the scenario ( fff.ij) . The black, red, and blue solid lines are for the 
components of <fe\, $2; and $3, respectively, and the black, red, and blue dashed lines for 
the components of§A u , ^A d , and§A e , respectively. 



CP- violating phases of the gaugino mass parameters (a = 1, 2, 3) and the dashed lines for 
those of the trilinear A parameters (a = 4, 5, 6). The dominant components are (^E dT1 ^ , 

(E dn ) , (E dH «J , fO dD ) , (O d A , and (o Acp j , reflecting strong dependences on 
the CP-violating phases of the gaugino mass parameters. The $1,2,3 components grow 
as tan/3 increases, implying that the EDM constraints on the CP phases $1,2,3 become 
stronger for larger tan/3. On the other hand, the &A u ,A d ,A e components are less than unity, 
except for O d °. This implies that the EDM constraints on these phases are weaker, and 
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Figure 4: The same as in Fig. but for expansions around the CP-conserving point <p a = 
180°. 

that the Deuteron EDM may be large enough to be observed in the proposed experiment. 
We also observe that the component (0 Acp ) > 10~ 2 when tan/3 > 30 and <p a = 0°, 
implying that it could give rise to the CP asymmetry larger than 1 % when $2,3 ~ 0° but 
$a„ is large, about 100°. 

In Fig. [51 we display the cosines C*o,E a of the relative angles between the three ob- 
servable vectors O dD ' dM,j4cp and the three EDM-constraint vectors E dT1 ,dn ,dH s. The upper 
frames are for the case of (p a = 0° and the lower ones for (p a = 180°. We see that the 
cosines between the observables and E dT1 (black solid lines) are reasonably small, except 
the case of the muon EDM. In the case, the Thallium EDM is dominated by the elec- 
tron EDM, resulting in the high degeneracy between the Thallium EDM constraint and 
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Figure 5: The cosines between the three observable vectors O dD,c! '" ylcp (left, middle, right 
frames) and the three EDM- constraint vectors E dT1 ' dn,dH s (solid, dashed, dash-dotted lines) 
as functions of tan /3 for ip a = 0° (upper) and ip a = 180° (lower). 

the muon EDM §. This alignment between the Thallium and muon EDMs in the scenario 
under consideration leads to a prediction of the muon EDM that is below the projected 
sensitivity, as we show in the next Section. On the other hand, the cosines between the 
Qd D ,Acp anc [ E dn (red dashed lines) are larger than those between O dD,j4cp and E dH s (blue 
dash-dotted lines) for (p a = 0°. In the case ip a = 180°, the cosines between the observable 
and the Mercury EDM-constraint vectors are slightly larger. 

Having the vectors representing the EDM constraints and observables in hand, one can 
combine them to construct the optimal directions in the 6D space of CP-violating phases, 



§The two EDMs are not exactly degenerate, due to the additional contribution to c?ti from the electron- 
Nucleon interaction Cs ei^e NN, which becomes larger as tan/3 increases. We find that the degeneracy 
is lifted by an amount of O(10 -4 ) for large tan/3, which it is too small to have visible effect in the figure. 
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using Eq. (12. 7ft . Specifically, we consider the three optimal directions that maximize dp>, 
dn, and Acp, respectively, taking into account the three existing EDM constraints. As 
comparisons to them, we also consider two other reference directions, which have A$x = 
= and A$ 2 = A$ 3 = 0, where A$ denotes the difference from the corresponding 
CP-conserving point. These two reference directions can be constructed by defining 

K = Me aM& ^E a p E\E^N^N^ , (3.6) 

where, for each direction, the two null directions Nj^' 2 ' are chosen as 

= (1, 0, 0, 0, 0, 0) , N® = (0, 0, 0, 0, 0, 1) for the direction A$ x = A$ Ae = 0, 
ArW = (o, 1, 0, 0, 0, 0) , N {2) = (0, 0, 1, 0, 0, 0) for the direction A$ 2 = A$ 3 = . 

(3.7) 

We show in Figs. |6] and [7] the absolute values of the six components of the five 
normalized vectors for (p a = 0° and (p a = 180°, respectively. We first observe that the 
$1,2,3 components (solid lines) are relatively small, and decrease as tan f3 increases. Hence, 
all the directions are mostly given by some combination of <&a u (black dashed line) and &A d 
(red dashed line) directions, with the exception of the A$2,3 = direction, which is mainly 
aligned with the $ J 4 e (blue dashed line) direction. The component is generally larger 
than the <&A d component in the optimal directions, except in the case of with tp a = 0° 
and tan /3 > 20, as seen in the middle-upper frame of Fig. [6j 

Finally, we consider the products $* • O of the 6D vectors of the normalized optimal 
directions and the observables. The products determine the sizes of the observables along 
the directions through the relations given in Eq. ( 12. 8 p when 0* = 1. As we see in the next 
Section, 0* could be as large as ~ 100 before the small-phase approximation breaks down 
and one of the three EDM constraints is violated. We show in Fig. [S] the products for the 
directions optimized for g? d , and Acp, which are denoted by the thick black solid, red 
dashed, and blue dotted lines, respectively. The thin lines are for the products along the 
reference A$i jJ 4 e = (thin black) and A$2,3 = (thin magenta) directions. The upper 
frames are for (p a = 0° and the lower ones for (p a = 180°. We observe that the directions 
constructed using the geometric prescription given in the previous Section do indeed give 
the optimal values of the observables, which are larger than those along the other optimal 
and reference directions, sometimes even much larger, depending on tan ft. 



4 Optimal Values of CP- Violating Observables in the 
MCPMFV SUSY Model 

As preparation for presenting our numerical results for g?d, d^, and Aqp, we first examine 
the magnitudes of the Thallium, neutron and Mercury EDMs along the three optimal and 
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<p I 

^ Ad 1 

A * 



Figure 6: The absolute values of the six components of the five normalized direction vectors 
considered in the text for (p a = 0°. The solid lines represent the $1,2,3 components, and the 
dashed lines the $A u ,A d ,A e components. 



two reference directions considered in the previous Section. As long as the small-phase 
approximation is valid, the EDMs should lie below the current experimental limits, but the 
approximation breaks down for large phases, leading to non-vanishing EDMs larger than 
the limits. The breakdown of the small-phase approximation would limit the maximum 
values of 0* in Eq. (12. 8 j) to values depending on the choice of direction and, accordingly, 
limit the maximum values of the CP-odd observables that can be found in this approach. 

Figure [9] shows the three constrained EDMs for tan ft = 40. From this Figure, one 
can read off the maximum allowed value of 0* for each direction. For example, from the 
upper frames with <p = 0°, (</>*) max ~ 70 , 50 , 70 , 50 , and 50 for the do-optimal, d M -optimal, 
Acp-optimal, A$i^ e = 0, and A$ 2 ,3 = directions, respectively, with the most important 
constraints being those provided by G?n g , <^H g , <feg, du g , and g?ti, respectively. Inserting the 
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values (0*) max into the corresponding products $* • O, one may get (near-)maximal values 
of the CP-odd observables in the linear regime. 

Fig. ITUl shows the three CP-odd observables along the five directions when tan /3 = 40. 
The experimental upper bounds on the Thallium, neutron and Mercury EDMs have been 
imposed. In the case of the Deuteron EDM with (p a = 0°, the linear regime ends at 0* ~ 
±30, and beyond this range it reaches saturation at a value ~ |9|, subsequently decreasing 
until the maximal value of the angle (0*)~^* oo ~ |70| is reached. When (p a = 180°, the 
linear regime ends at <p* ~ ±30. Beyond this point, the Deuteron EDM along the g?d 
optimal direction saturates, but it continues to increase up to ~ |6| along the and 
A^i,^,, = directions. These two examples demonstrate that the existing EDM constraints 
do not exclude the observation of d^. 

On the other hand, the muon EDM is always below the projected sensitivity, since 
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: d D dir 



A *i,Ae =0 dir - — : A *2,3=° dir - 




Figure 8: The products $* • O along the directions optimized for djy, d^ and Aq P , which 
are denoted by the thick black solid, red dashed and blue dotted lines, respectively. The thin 
lines are for the products along the A$ 1]J 4 e = (thin black) and A$2,3 = (thin magenta) 
directions. The upper frames are for tp a = 0° and the lower ones for (p a = 180°. 



$* ■ O dt± is forced to be small by the (near-) degeneracy with the Thallium EDM, see Fig. |HJ 

The CP asymmetry Acv{b — > 57) has a larger linear regime, extending over the almost 
whole range with |0*| < 100, and could as large as ~ 2% when ip = 0°. Such a value could 
be observed at a future Super B Factory [21]. Unfortunately, the case with tan /3 ~ 40 gives 
too small a branching ratio B(b — >■ sj) ~ 1 x 10~ 4 B In order to respect the constraint 
from B{b — > S7), we consider the case of small tan/3 = 10 with ip = 0°, see Fig. [TTJ This 
case is compatible with the 5(6—7- sj) constraint at the 2-cx level, but the predicted A CP 
is too small to be observed. On the other hand, the attainable values of g?d and <i M are 

^See also [T2], in which the same scenario has been analyzed. 
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Figure 9: The values of the three EDMs along the three optimal and two reference directions 
for tan/3 = 40, with (p = 0° (upper) and (p = 180° (upper). The line styles are the same as 
in Fig. [21 with the additional horizontal lines representing the EDM constraints. 



similar to those in the case when tan (3 = 40. 

Finally, we show in Fig. [12] the 6 CP- violating phases at the GUT scale and the 3 
CP-violating phases of the third-generation A parameters at the SUSY scale when (p a = 0° 
and tan/3 = 10. We observe that the CP-odd gaugino phases $i, $2 and $3 can only be 
as large as 4°, 1° and 2°, respectively, whereas the CP-odd trilinear A-phases &A u ,A d ,A e at 
the GUT scale can be as large as ±70°. These CP- violating phases are suppressed at the 
SUSY scale by RG running from the GUT scale [12] , but sizeable non-trivial CP- violating 
phases are still allowed at the SUSY scale: A$ At ~ ±4°, A$ Ab ~ ±10°, and A$ At ~ ±40°. 
Similar magnitudes of the CP- violating phases are attainable for (p a = 180°. 
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Figure 10: The three CP-odd observables along the three optimal and two reference direc- 
tions for tan/3 = 40, with (p = 0° (upper) and (p = 180° (upper). The line styles are the 
same as in Fig. 0. 

5 The 7D Extension Including #qcd 

In the previous Sections, we implicitly assumed that the CP-violating QCD #-term: 



CY — ~ 

8tt ^ 



(5.1; 



vanishes, where G' iu,a = e tJ-upcr Gp cr /2 and the parameter 9 is given by the sum of the QCD 
#qcd and the strong chiral phase for the quark mass matrix as 



9 = #qcd + Arg Det M q . 



(5.2) 



In the weak basis where Arg Det M q = 0, we have 9 = #qcd- The QCD #-term (15. ip would 
be set to zero, e.g., if there is a QCD axion in the theory. Otherwise, the dimension-four 
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Figure 11: The same as in Fig. [Wj but for tan /3 = 10 and (p = 0°. 



operator ( 15. ip would in general contribute to the neutron, Mercury and Deuteron EDMs, 
e.g., through the CP-odd pion-nucleon-nucleon interactions 

C wNN = g^ N Nr a N n a + £ ] NN ~NN ix° . (5.3) 

QCD sum rule techniques have been used to estimate the contribution of the 9 term to the 
neutron EDM in [22l[23] ft 



d n {9) = (0.4 ± 0.2) [x(4e d - e u )mj 

= ^ 1±0 - 5 )(2ileW' X2 - 5><10 " 16eCm ' (5 - 4) 

where the reduced mass m* = m u md/ {m u + mA = — (1/3) e, e u = (2/3) e, and the 
condensate susceptibility x = — 5.7±0.6 GeV~ 2 o The d n {9) contribution is to be added to 
the other EDM contributions from the electric and chromoelectric dipole moments of light 
quarks, the Weinberg operators, and the CP-odd four-fermion interactions that are induced 
by the CP- violating phases in the soft SUSY-breaking sector, i.e., by the six independent 
MCPMFV CP phases considered above [23]. 

The leading contribution of the 9 term to the Mercury EDM is expected to be through 



"Though there are some inconsistencies between the two references, the second numerical equation 
in [23] is not affected by them. 

** It has been argued that, in the presence of an axion, 9 should be replaced by 9- m d, which is determined 
dynamically via the chromoelectric dipole moments of the up, down, and strange quarks d 1 ^ d s . However, 
since we lack full knowledge of all the relevant higher-order corrections [H] , we will neglect 6*i n d contributions 
to EDMs in our analysis. 
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Figure 12: The 6 CP-violating phases at the GUT scale and the 3 CP-violating phases of 
the third- generation A parameters at the SUSY scale along the 3 EDM-free directions and 
2 reference directions when tp a = 0° and tan/3 = 10. 

its contribution to the pion-nucleon-nucleon interaction coefficient <?£jviv : 

d Hg (0) = +(1-8 x 1(T 3 GeV" 1 ) egV N (9) . (5.5) 

The 9 contribution to the coupling is suppressed by the factor (ma — m u )/m s , and given 
by (HI 

in _ nm f) m 7 — m _ _ _ 

(5.6) 



g^ N {9) = ^ m \ mu (N\uu + dd - 2ss\N) ~ 1.1 x 10~ 3 9 



which results in 



d Hg {9) ~ +2.0 x lO^eGeV- 1 ~ 3.9 x 10~ 20 ^ecm. 



(5.7) 
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The leading contributions of the 9 term to the Deuteron EDM are given by [TT] 

d D {6) ~ -e [(3.5 ± 1.4) + (1.4 ± 0.4)] x 10~ 3 9 GeV" 1 ~ -9.7 x 10~ 17 9ecm. (5.8) 

The first term is the leading-order QCD sum-rule estimate of the 9 contribution to the 
sum of the proton and neutron EDMs, which enters only via subleading isospin-violating 
corrections. The second term arises from the coupling g^ N (9): 

d«™(9) = - e9 ^£^f ] 7^|i - -(LSiO^e^GeV- 1 , (5.9) 

where g n NN — 13.45 and £ = y/m p e/m n , with e = 2.23 MeV being the Deuteron binding 
energy. 

Including non-vanishing 9 together with the six MCPMFV CP phases $1,2,3 and 
*&A u ,A d ,A e i there is a total of seven CP phases. The six-dimensional geometric construction 
of optimal EDM-free directions in Section 2 can easily be extended by one more dimension 
in this case. As in the 6D case, the 7D EDM-constraint and observable vectors are given 
by E = VE and O = VO, but with 



d d d d d d d 



<9V <9$ 2 ' <9$ 3 d$ A J d$ A ' d$ A ; 39 



(5.10) 



As before the CP-violating phases $1,2,3 and §A u ,A d ,A e are specified in degrees and we 
normalize 9 in units of 10~ 10 : 

9 = 9xl0 10 . (5.11) 

With this normalization, with 9 = 1, we have d n (9) = 2.5 x 10 _26 ecm which is very near 
to the current experimental bound d^ xp = 3 x 10~ 26 ecm. 

In Fig. [131 we show the absolute values of the components of the three 7D EDM- 
constraint and the three 7D observable vectors in the linearized CP-violating version of the 
scenario ( 13. ip . taking <y5 a= i_ 6 = 0°. Comparing with Fig. [3]of the 6D case, we see that: (i) 
E dT1 , O d ", and Acp are unchanged, with vanishing seventh 9 components, (ii) the first 
six components of E dn , E dfi s and O dn are unchanged, and we have the new 9 components, 
(E dn ) 7 ~ 0.83, (E dH °) 7 ~ 0.13 and (e 6 ^ ~ -3.2, which are independent of tan (3. 

Having obtained the 7D E and O vectors, one can construct the optimal EDM-free 
direction maximizing O as in the 6D case: 

K = M e aM s»u P E d ^ Es Hs B^ vp , (5.12) 

where the 3-form B^ vp is given by 

B^p = e^pxvru, O x E d a ^ E d - E*** . (5. 13) 



22 




tanP tanP tanP 




tanP tanP tanP 

Figure 13: The absolute values of the components of the three ID EDM- constraint vectors 
(upper) and those of the three ID observable vectors (lower) in linear expansions around 
(p a= i_Q = 0° , as functions of tan /3 for the scenario $3. The magenta dashed lines 
represent the 9 = 9 x 10 10 component, while the other lines are the same as in Fig. 

We show in Fig. [TH the absolute values of the seven components of the three optimum 
directions (upper) and the two reference directions with A$! = A$^ e =9 = (lower left) 
and A$2 = A$3 = 9 = (lower middle), considering the scenario (13. ip with (f a =i-6 = 0°. 
Comparing with Fig. [6] of the 6D case, we see that: (i) the vectors in the two reference 
directions with 9 = remain the same, (ii) the dj)- and ci^-optimal directions can now have 
sizeable $ 12 components, while the 9 component dominates when tan/3 > 10, (Hi) the 
Acp-optimal direction is still dominated by the $a„ component when tan /3 > 7. Fig. [T5l 
shows the products with the three observable vectors of the three optimal EDM-free direc- 
tions and the two reference directions for the same scenario. We observe again that the 
optimal direction found using our geometric construction gives the largest value for each 
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Figure 14: The absolute values of the seven components of the five normalized direction 
vectors for the 7D CP-violating version of the scenario 113. in a linear expansion around 
<y5 Q=1 _ 6 = 0°. The lines are the same as in Fig. 0, with additional magenta lines for the 9 
components. 



corresponding observable. Comparing with Fig. [8] (upper frames) of the 6D case, the prod- 
ucts can be larger by more than an order of magnitude for the Deuteron and muon EDMs. 
On the other hand, they remain more or less the same for Acp, due to the dominance of 
the $a„ component in the ^cp-optimal direction. 

In the upper frames of Fig. [TBI we show the three EDM constraints, assuming tan /3 = 
10 and <^ Q= i_6 = 0°. As in the case of Fig. [§1 one can read off the maximum values of 
<fr* for each EDM-free direction from the figure: (0*) max ~ 25 , 25 , 50 , 40 , and 45 for the 
do-optimal, ^-optimal, Acp-optimal, A$i t A e = 9 = 0, and A$ 2 3 = 9 = directions, 
which are mainly constrained by dxi, ^ti, ^ti, ^H g , and d-pi, respectively. Multiplying the 
values (</>*) max to the corresponding products $* ■ O shown in Fig. (TSJ we find the maximum 
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Figure 15: The products $* ■ O along the three dp,-, d^-, and A^-p- optimal and the two 
arbitrary directions for the ID CP-violating version of the scenario Ii3. 1\) with y5 Q= i_6 = 0°. 
The lines are the same as in Fig. 

values of the CP-odd observables in the linear regime, which are shown in the lower frames 
of Fig. [16] Comparing with Fig. [11] of the 6D case, we see that the maximal value of 
the Deuteron EDM is greatly enhanced, becoming as large as ~ 70 times the projected 
sensitivity. On the other hand, the muon EDM is still below the projected sensitivity, and 
the the maximal value of the CP asymmetry Acp(b — > sj) could be only ~ 0.08%. 

Finally, in Fig. [T71 we show the 6 CP-violating phases at the GUT scale and the 3 
CP phases of the third-generation A parameters at the SUSY scale along the 3 EDM-free 
and the 2 reference directions in the 7D case. Comparing with Fig. [12] for the 6D case, 
we observe in the top panels that $i and $2 can be substantially larger, as large as ~ 20° 
and ~ 5°, respectively, for (0*) mdx ~ 25 along the dp,- and <i M -optimal directions denoted 
by the thick solid and dashed lines. We also note in the middle and bottom panels that 
the phases of Ad >u ,e are somewhat larger than in the 6D case. Finally, we note (not shown) 
that 9 could be as large as ~ 5 x 10 -9 along the Acp-optimal direction. 
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Figure 16: Upper: The three EDM constraints along the three EDM-free and two reference 
directions for the scenario A3. 1\) with tan/3 = 10 and (p a =i-e = 0°. Lower: The three 
CP-odd observables in the same scenario. The lines are the same as in Figs. andUU 



6 Conclusions 



We have proposed in this paper a novel geometric technique for optimizing the possible 
values of CP-violating observables in the presence of the strong constraints due to upper 
limits on EDMs. Our geometric approach enables us to separate the EDM-free subspace off 
the full CP-phase parameter space in the linear approximation. Knowing the parametric 
dependence of a given observable, we can analytically construct the extremal direction in 
the full parameter space, along which the observable gets maximized. Since our approach 
is analytic, it becomes exact in the linear approximation and is much more efficient and 
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Figure 17: TTie # CP phases at the GUT scale and the 3 CP-violating phases of the third- 
generation A parameters at the SUSY scale along the 3 EDM-free and 2 reference directions. 
The scenario 113. 1\) is taken with tan/3 = 10 and y5 a= i_6 — 0°. The line styles are the same 
as in Fig. [XT 



accurate than the naive type of scan that is usually made in literature [T7] |H| 

We have demonstrated the applicability of this technique in two cases: the 6D case 
of the MCPMFV version of the MSSM, and a 7D extension to include the QCD vacuum 
phase. We have illustrated this approach within a class of CP-violating models that extend 
and generalize the well-studied SPSla benchmark point of the CMSSM. For any specific 
benchmark point, the values of CP- violating observables are in general bounded in magni- 

^For comparison, we note that, for a theory with 6 free parameters, a search within a 6-dimensional grid 
with 100 points in each coordinate would require 10 12 scan points. Instead, our geometric method only 
involves straightforward sums over a 6-dimensional Levi-Civita tensor with 6 ! = 720 non-zero components, 
along a radial line of 50 points. 
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tude, since the ranges of the CP- violating phases fa are all compact: fa G [0, 2ir). Based on 
linear expansions of CP- violating observables around CP-conserving points with fa = 0, n, 
our approach gives rather accurate estimates of the true maximal values of the CP- violating 
observables. 

Using this approach in the 6D MCPMFV case, we find values of the Deuteron EDM 
that may be an order of magnitude larger than the prospective experimental sensitivity. 
This range is increased by almost another order of magnitude if our optimal geometric 
construction is extended to the 7D case that includes the CP-violating QCD vacuum phase. 
Hence, the Deuteron EDM may provide indirect useful information about the possible 
presence of a non-vanishing QCD vacuum phase, complementing the current experiments 
on EDMs and the experimental searches for axions [26] . On the other hand, we find that the 
maximal values of the muon EDM are somewhat below the likely experimental sensitivity 
in both the scenarios with and without the QCD phase. We also find that the CP-violating 
b — > S7 decay asymmetry Aqp is too small to be observed, once the stringent constraint 
from B(b — >■ 57) is taken into account. Likewise, the I^-mixing phase <pB a turns out to be 
close to the small SM value in both the scenarios studied. 

Our geometric approach could easily be extended to other supersymmetric scenarios. 
For example, we have not made a systematic survey of all the possibilities in the MCPMFV 
SUSY model that arise as generalizations of other CP-conserving benchmarks. It may also 
be interesting to extend this approach to a wider class of CP-violating models within the 
general MSSM framework. Since it has many more CP-violating parameters, our approach 
may be a useful guide to the possibilities opened up by this larger parameter space. Needless 
to say, our approach may also be interesting for other scenarios for CP violation, beyond 
the MSSM and indeed supersymmetry. Finally, we observe that this geometric approach is 
not restricted to issues related only to CP violation. It could find broader applicability to 
other problems where one wants to maximize observables subject to a set of constraints. 
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A Two-Loop Gaugino Contributions to EDMs 



In this Appendix we calculate a particular set of EDM contributions induced by two-loop 
diagrams involving chargino (x*) and neutralino quantum effects. The relevant tree- 
level interactions are: 

C w ± m = (9w + ^jP L + 9w + ^j?n) XJ + H.c, (A.f) 

where the f/^-boson couplings are given in the CPsuperH manual [27], and the W /± -boson 
couplings are given by 

gK +mj = -N: 4 (C R )* 2 + V2N; 2 (C R y n , (A.2) 

in terms of chargino and neutralino mixing matrices. The charged-Higgs couplings to the 
SM particles are given by 



C H ± ff , = - 9ff ,H + f [g H +f f , +i%g H +ff>) f + H - c -> ( A - 3 ) 

where 

Qui = > 9 S H+ m = t fi /2 , g p H+9l = -it p /2 . (A.4) 
We include the threshold corrections due to the exchanges of gluinos: 
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(A.5) 



where 

A, = -^^M;i{M% M|, |M 3 | 2 ) . (A.6) 

We note that this threshold correction induces a dependence on the gluino mass phase $3. 
At the two-loop level, charginos and neutralinos induce CP-violating interactions in the 
7-if ± -W /=F and ■y-W ± -W T couplings, which in turn produce non-zero electron and quark 
EDMs. 
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A.l EDM Induced by the j-H^W* Coupling 

In addition to squarks [28], charginos and neutralinos may also induce a complex and CP- 
violating effective , y-H ± -W T coupling [19]. The CP- violating effective r y-H ± -W T coupling 
can then give rise to electron and <i-quark EDMs. Specifically, defining as f = l,d all 
fermions with weak isospin T z = —1/2, we have 
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(A.7) 



AB 



'w+x^xj y W+x iX~ 



with A, B = L, R. The loop function J (a, b) is defined in terms of the function J(x), 

, v *jC In oc , . _ s 

J x = , A.9 

x — 1 

by J (a, b) = [J (a) — J(b)]f (a — b). For a = b, the loop function J (a, b) takes on the simple 
form: J(a, a) = (— In a + a — l)/(a — l) 2 , with J(l, 1) = 1/2. Observe that the expression 
9 H+jf + ^9 H+ff = (t/3 m d/'m'u) for /' = Z (d) is real at the tree level, but becomes in 
general complex when gluino threshold corrections are included, as discussed above. 



A.2 EDM Induced by the j-W^W^ Coupling 

Quantum loops of charginos and neutralinos generate P- and CP-odd interactions in the 
, y-W ± -W T coupling [29H3T] , which in turn produce non-zero fermion EDMs at the two-loop 
level. The analytic results we use are based on the latest calculation in [21]. In detail, if 
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we define / = l,d, the contribution of the , y-W ± -W T coupling to the /-particle EDM is 
given by 



'dJ_\ WW 



a 

Xsm. 



J J Wr 4 , tv) , (A.10) 



where 



with rj = m?±/M^ and = m~o/M^, and J(a, 6) being defined after (|A.9[) 



A. 3 CPsuperH2 . 2 Interface 

The additional two-loop EDMs discussed above have been implemented in the public version 
of the CPsuperH2 . 2 code. The output of the two-loop EDM calculations and the strange- 
quark chromo-EDM is contained in the auxiliary array RAUX_H. The additional two-loop 
EDMs involving charginos and neutralinos have been added to the corresponding EDMs 
presented in Ref. [25] to yield the total EDMs. Specifically, the following assignment of 
variables (all in units of cm) has been made in the updated code CPsuperH2 . 2: 

• Two-loop electron EDMs induced by the r )-H ± -W Zf and r )-W ± -W T couplings: 

RAUX_H(205) = (dfje) wm , RAUX_H(206) = (df /e) ww . (A.12) 

• Two-loop down-quark EDMs induced by the ^-H^-W^ and ^-W^-W 11 couplings: 

RAUX_H(225) = (df /e) WH * , RAUX_H(226) = (df/e) WW . (A.13) 

• The corresponding two-loop strange-quark EDMs: 

RAUX_H(235) = (df /e) wm , RAUX_H(236) = (df /e) ww . (A. 14) 

• The chromo-EDM of the s-quark: 

RAUX_H(400) = d c s = (df )** + (df)*° + (df) 9 + (d c s ) H , (A. 15) 

where the individual contributions are given by 

RAUX_H(401) = (d c s )** , RAUX_H(402) = (df )*° , 

RAUX_H(403) = (dff, RAUX_H(404) = (d c s ) H . (A.16) 
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